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O i ABSTRACT 

O 

■ About 30-40 percent of the expected number of baryons is still missing in 

the local Universe (z<0A). They are predicted to be hiding in a tenuous web 
of intergalactic gas at temperatures of about 10 5 — 10 7 K (the Warm-Hot In- 
U- tergalactic Medium: WHIM), direct consequence of the collapse of matter 

i-C ■ into dark-matter potential wells during the continuous process of structure for- 

mation. Detecting this matter has had limited success so far, because of its 
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low-density and high temperature, which makes it difficult to detect with current 
far-ultraviolet and X-ray instrumentation. Large observational efforts and fine- 
tuned observing strategies are needed to finally dig out these hot baryons, decide 
on the existence of this matter in a WHIM and, by doing so, test fundamental 
theoretical predictions. 

Here we present the first results from our pilot 500 ks Chandra-LETG Large 
Program observation of the soft X-ray brightest source in the z > 0.4 sky, the blazar 
1ES 1553+113. We identify a total of 11 possible absorption lines, with single- 
line statistical significances between 2.2-4.1(7. Six of these lines are detected at 
high single-line statistical significance (3.6 < a < 4.1), while the remaining five 
are regarded as marginal detections in association with either other X-ray lines 
detected at higher significance and/or Far-Ultraviolet (FUV) signposts. Three of 
these lines are consistent with metal absorption at z ~ 0, and we identify them 
with Galactic 01 and CII. The remaining 8 lines may be imprinted by inter- 
vening absorbers and are all consistent with being high-ionization counterparts 
of Far-Ultraviolet (FUV) HI and/or OVI IGM signposts. In particular, four of 
these eight possible intervening absoption lines (4.1a, 4.1er, 3.8<r and 2.7a), are 
identified as CV and CVI Ka absorbers belonging to two WHIM systems at 
zx = 0.312 and (zx) = 0.133, which also produce broad HI and OVI absorption 
in the FUV. The true statistical significances of these two X-ray absorption sys- 
tems, after properly accounting for the number of redshift trials, are 5.8a and 
3.8a, respectively. 

Subject headings: WHIM, Absorption Lines 

1. Introduction 

According to the latest baryon census in the local (z<0A) Universe (Shull, Smith & Dan- 
forth, 2012 and references therein), about 30-40% of the total number of baryons predicted 
by Big-Bang Nucleosynthesis (e.g. Kirkman et al. 2003), inferred by density fluctuations of 
the Cosmic Microwawe Background (e.g. Bennett et al. 2003; Spergel et al. 2007), and seen 
at z ~ 3 in the "Lya Forest" (e.g. Rauch, 1998; Weinberg et al. 1997), are missing at the 
present epoch. 

Of the \ow-z baryons we can account for in a recent census (Shull, Smith & Danforth, 
2012), (18 ± 4)% lie in a collapsed phase (stars, cold atomic and molecular gas in galaxies, 
circumgalactic medium - CGM -, and hot virialized gas in groups and clusters of galaxies), 
(28 ± 11)% lie in the local photoionized Lya forest, and finally, as much as (25 ± 8)% , are 
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estimated to be found in the OVI/BLA-traced Warm-Hot Intergalactic Medium Q| . The 
remaining (29 ± 13)% of baryons still remains elusive, and are likely to lie in the hot phase 
of the WHIM, at temperatures larger than logT~ 5.5 K. 

This is predicted by nearly all hydrodynamical simulations for the formation of struc- 
tures in a A-CDM Universe (e.g. Cen & Ostriker, 2006, Tornatore et al. 2010, Bertone et al. 
2010, Branchini et al. 2009, Smith et al. 2011). According to simulations, starting at red- 
shift of z ~ 2 some of the IGM was shock-heated to temperatures of 10 5 — 10 7 K during the 
continued process of collapse and structure formation, and enriched up to Z WHIM = 0.1 — 1 
Z by galaxy super-winds (GSW, e.g. Cen & Ostriker, 2006) already at redshifts of ~ 0.5, 
remaining essentially unchanged from then to z ~ 0. The first pieces of observational evi- 
dence of this violent process of structure formation are starting to be gathered in the densest 
nodes of this filamentary IGM, at and inside the virial radius of massive clusters of galaxies 
(e.g. Keshet et al. 2012). However, diffuse WHIM filaments are too tenuous to be detected 
through their bremsstrahlung and line emission with current Far- Ultraviolet (FUV) or X-ray 
instruments (e.g. Yoshikawa et al. 2003). Absorption is more promising. 

The lowest temperature tail of the WHIM distribution, from T~ 10 5 up to ~ 10 5 5 K 
(where the fraction of OVI peaks), has been most likely detected in the FUV, with FUSE and 
the STIS spectrometer onoard the Hubble Space Telescope (HST), in OVI and narrow Lyct 
absorbers, though the fraction of this gas that actually belongs to the shock-heated WHIM 
is highly uncertain (e.g. Danforth & Shull, 2008). More recently, the advent of the larger 
throughput Cosmic Origins Spectrograph of HST (COS, hereinafter) has enabled the detec- 
tion of somewhat hotter (T~ 10 5 — 10 5 ' 7 ) and possibly truly shock-heated WHIM, through the 
detectability of the strongest (i.e. large column density) Broad Lyman a absorbers (BLAs) 
with thermal Doppler parameter b t h = y / 2kT/m p ~ 40 — 90 km s^ 1 (Danforth, Stocke & 
Shull, 2010) § . These BLAs are in principle good tracers of the warm (T~ 10 5 — 10 5 ' 7 K) 
portion of the WHIM, where a significant fraction of the WHIM (~ 30 — 40%, e.g. Cen & 
Ostriker, 2006) is supposed to reside. However, observing b t h — 40 — 90 km s _1 BLAs re- 



1 This estimate is based on an average ionization and metallicity, versus column density, correction derived 
by hydrodynamical simulations, and therefore suffers not only the large scatter found in the given simulation 
between these parameters (i.e. Fig. 8 in Shull, Smith & Danforth, 2012), but also the scatter of theoretical 
predictions from simulation to simulations (e.g. Cen & Ostriker, 2006, Tornatore et al. 2010, Bertone et al. 
2010, Branchini et al. 2009, Smith et al. 2011). The resulting number of baryons still missing also depends 
on these uncertainties. 

2 The thermal Doppler parameter, is typically narrower than the measured Doppler parameter 6, which 
includes a combination of internal turbulence and Hubble-flow broadening. Danforth, Stocke & Shull, 2010, 
estimate an average correction factor (b/b t h) — 1-2. 
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quires data of excellent quality and accurate continuum definition to detect the low-contrast 
features. In addition, measuring HI without an ionization correction [f| will always underesti- 
mate the actual baryon fraction in the fraction of hot IGM sampled by these tracers. Given 
the temperature of this gas, such an ionization correction can only come from the X-ray 
detection of the highly ionized counterparts to such BLAs. Finally, the T>10 5 ' 7 K portion 
of the WHIM (the vast majority of the WHIM mass) can only be detected through highly 
ionized metals in the X-rays. 

However, detecting the metal tracers of the bulk of the "Missing Baryons" in the "X-Ray 
Forest", has proven to be extremely difficult. This is because of the unfortunate combination 
of (a) the still limited resolution (R ~ 400 at the location of the OVII Ka transition, at 
A = 21.6 A) and the low throughput (A-^g ~ 15 — 70 cm 2 ) of the current high-resolution X-ray 
spectrometers [the XMM-Newton Reflection Grating Spectrometer (RGS, den Herder, et al. 
2001) and the Chandra Low Energy Transmission Grating (LETG: Brinkman et al. 2000)]; 

(b) the lack of bright (/o.5-2fceV > 10 -11 erg s -1 cm -2 ) extragalactic point-like targets; and, 

(c) the dramatic steepening in slope (Aa> 1.5) of the predicted number density of metal He- 
like WHIM filaments per unit redshift at ion column densities > 10 15 cm -2 . Current evidence 
is still limited, and highly controversial. 

The statistically strongest evidence to date comes from the evidence for large amounts 
of baryonic matter in two intervening WHIM filaments along the line of sight to the nearby 
(z = 0.03) blazar Mkn 421 (Nicastro et al. 2005a,b). However, these results are controversial 
(Kaastra et al. 2006, Rasmussen et al. 2007, Danforth et al. 2011, Yao et al. 2012; but see 
also Williams et al. 2006, Nicastro et al. 2008, Williams et al. 2010). 

Over the last few years there have been new discoveries reported (Buote et al. 2009: 
B09; Fang et al. 2010: F10; Nicastro et al. 2010, Zappacosta et al. 2010, Zappacosta et al. 
2012). However, such tentative discoveries are either serendipitous (and so biased toward 
the strongest systems: Nicastro et al. 2010; Zappacosta et al, 2012), or trace extreme galaxy 
overdensity regions (Zappacosta et al. 2010) or even just galaxy halos (Buote et al. 2009; 
Zappacosta et al. 2010; Williams et al. 2012). In all cases, these sample only the hottest 
and densest IGM or circumgalactic medium, which is not representative of the majority of 
the WHIM. 

Thus a firm, high-significance detection of the bulk of the missing baryon mass has 



3 The measured HI Doppler parameter gives only an upper limit on the gas temperrature, and in addition, 
the proper ionization correction depends also on the gas density because of the second-order contribution of 
photo-ionization, and this can only be estimated through multi-ion metal measurements, e.g. Nicastro et al. 
2002). 



-5 - 



not been secured yet. To do so, the only pursuable observational strategy, with current 
instrumentation, is to first identify "true" gaseous signposts (e.g. Mathur et al. 2003) of 
the T~ 10 5 3 — 10 5 ' 7 WHIM (i.e. BLAs), and then search for their metal counterparts in 
the soft X-rays, along optimally selected sightlines. Over the past two years we selected 
1ES 1553+113 as "the best WHIM target in the Universe' (because of its extreme soft 
X-ray brightness, relatively high redshift and the presence of BLA signposts in its HST- 
COS spectrum - Danforth et al. 2010), and started the X-ray campaign with a pilot 500 ks 
Chandra cycle 12 observation of 1ES 1553+113. In our observing proposal, we conservatively 
tuned the exposure of this first pilot observation to test the reliability of the strongest FUV 
systems as gaseous signposts for X-ray WHIM filaments. The results are well beyond the 
most optimistic expectations and confirm the effectiveness of such an observational strategy. 
In this article we report on these results from the pilot Chandra-LETG spectrum. 

The paper is organized as follow. In §2 we describe the target of our LETG observation, 
and summarize the relevant results from the HST-COS observations reported by Danforth et 
al. (2010). Section 3 describes the Chandra-LETG data reduction and analysis, while in §4, 
we discuss our findings and present possible absorption system identifications. Finally in §5, 
we summarize our conclusions. In a forthcoming article ("paper II"), we will present a self- 
consistent modeling of the spectral features detected in the LETG spectrum of 1ES 1553+113 
by making use of our WHIM hybrid collisional-ionization plus photoionization models, and 
so deriving preliminary estimates of metallicities of these systems, the ionization-corrected 
baryon cosmological mass density of the WHIM at T~ 10 5 5 K, the number density of CV 
and OVII WHIM absorbers as a function of the ion column density. 



2. 1ES 1553+113 and the WHIM Signposts along its Line of Sight 

Four years ago, we selected 1ES 1553+113 as the optimal target for FUV and X-ray 
WHIM studies, because of its relatively high estimated redshift prior to spectroscopic de- 
terminations (z>0.25, from the non-detection of the host galaxy in a deep HST R image: 
Treves et al. 2007), its high Galactic latitude (b ~ 44°) and its extreme average 0.1-2.4 keV 
brightness of 2 mCrab (4 x 10 -11 erg s _1 cur 2 ), coupled with a bright and stable FUV flux 
of vf u (1600A) ~ 3 x 10" 11 erg s" 1 cm" 2 . 

1ES 1553+113 was observed with the HST-COS on September 22, 2009, for 3 orbits 
(1 orbit with the G130M and 2 with the G160M grating: Danforth et al. 2010). Here we 
summarize the main findings from this HST-COS observation, relevant to this work. The 
HST-COS spectrum of 1ES 1553+113 secured a robust lower limit on the spectroscopic 
redshift of 1ES 1553+113: Narrow HI Lja (NLA) and OVI absorption is observed at a 
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combined 19.3<r significance at z — 0.395 (Fig. 2 of Danforth et al. 2010), so establishing 
a spectroscopic lower limit on the redshift of 1ES 1553+113 of z>0A (almost twice the 
non-spectroscopic value upon which we based our original selection). Moreover, Danforth 
et al. (2010) report on the presence of 5 possible BLAs, at z = 0.04281, 0.10230, 0.12325, 
0.13334 and 0.15234, plus a rare and complex triple-HI-metal absorption system at z = 
0.1864 — 0.1899, spanning a range of 1200 km s -1 in velocity (Fig. 2, Table 1, in Danforth 
et al. 2010). We use these systems as gaseous signposts for either WHIM filaments (BLAs), 
or highly ionized galaxy outflows and hot circumgalactic medium from nearby galaxies (the 
triple-HI-metal system). 

Over the past 13 years, 1ES 1553+113 has been observed with an average 0.1-2.4 keV 
flux of 2 mCrab (1 mCrab = 2xl0~ n erg s" 1 cm" 2 ). A recent short (22 ks) XMM DDT 
and our 3-visit 500 ks LETG observation, confirmed this historical average, all recording 
Fo.i~2.4keV — 2 mCrab. This, combined with the redshift lower limit established by the COS 
spectrum for 1ES 1553+113, makes this target by far brightest, non-transient, X-ray source 
in the z > 0.4 sky. 



3. The Chandra-LETG spectrum of 1ES 1553+113 

Chandra observed 1ES 1553+113 between May- June 2011, for 500 ks, with the LETG 
spectrometer (constant resolution AA = 50 mA, i.e. resolving power R= 400 at A = 20 A) 
dispersed across the High Resolution Camera for Spectroscopy (HRC-S, Murray et al. 1985) 
plates. The observation consisted of three visits with exposures of 166.3 ks (2011, May 4-6), 
175.4 ks (2011, May 6-9) and 153.9 ks (2011, June 18-20). 



3.1. Data Reduction 

The data were reduced with the latest version of the Chandra Interactive Analysis of 
Observation software (CIAO v. 4.4.1, CALDB v. 4.5.2), following the standard processing 
procedures outlined in the HRC-S-LETG grating analysis thread □ . 

The Chandra HRC-S detector has virtually no spectral resolution, which prevents the 
separation of the different spectral orders dispersed by the LETG. For each observation, we 
therefore extracted all-order negative and positive source and background spectra with the 
CIAO tool tg_extract. This produced a total of six LETG spectra, three (one per LETG 



4 http://cxc. harvard.edu/ciao/guides/gspec_hrcsletg. html 
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visit) negative-order and three positive-order source and background spectra, extracted at a 
binning pace of AA = 12.5 mA, i.e. oversampling the LETG FWHM resolution of 50 mA by 
a factor of 4. 

Due to the impossibility of separating HRC-LETG spectral orders, the spectral mod- 
eling of these spectra must be performed by folding the fitting models with the sum of 
the convolution products of the redistribution matrices (RMFs) and ancillary response files 
(ARFs) belonging to the first N positive and negative orders, respectively. For our LETG 
observations of 1ES 1553+113, we conservatively decided to build our final all-order response 
matrix by adding up positive and negative orders up to N = 10. For each order, we used the 
CIAO tools mkgarf and mkgrmf 'to build 20 effective area and photon-redistribution matrices, 
which were then read together with the data into the CIAO fitting-package Sherpa, through 
the commands load_multi_arfs and load_multi_rmfs. 

To maximize the S/N of our spectra (we note that the HRC-S-LETG combination has 
an effective area of only A E fr — 15 cm 2 at both A = 22 and 44 A), we first co-added 
the three negative and three positive-order spectra (source and background) from the three 
observations (CIAO tool add_grating spectra), to produce 2 LETG spectra (one negative- and 
one positive-order spectrum), with net exposures of 495.6 ks each. Effective area matrices 
of the three observations, for each of the 10 negative and ten positive orders, were averaged 
together, producing 20 averaged effective area matrices (ARFs). Finally, we used the CIAO 
tool add-grating_orders to co-add the two negative and positive order source and background 
spectra and ARFs, to end up with a single LETG spectrum, containing all counts from the 
ten negative-order and the ten positive-order spectra of the three Chandra observations of 
1ES 1553+113. The final co-added spectrum, has roughly 250-300 net (i.e. background 
subtracted) counts per AA = 50 mA, resolution element, over the two interesting spectral 
bands covering the z < 0.4 OVII Ka (A = 21.6 - 30 A) and CV Ka (A = 40.2 - 56 A) 
ranges at redshifts where WHIM absorption might be found, for a signal-to-noise ratio per 
resolution element (S/N)~ 16. 

Because of our co-adding procedure, particular care must be taken to account for the 
wavelength calibration inaccuracy. The Chandra LETG-HRC-S spectrometer suffers a vari- 
able intrinsic wavelength uncertainty owing to the non-linearity in the effective dispersion 
relation. This non-linearity is due to a spatial imaging non-linearity of the HRC-S detector 
caused by the combination of (1) partial charge loss in the 3-tap read-out algorithm of each 
of the three HRC-S microchannel plates, each of which is divided into 190 taps along the 
dispersion direction, and (2) a hardware problem that corrupts the data from the position 
taps under a specific set of conditions. The net result is that photon event positions obtained 
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from processing HRC-S telemetry can be offset from their true values § . The wavelength 
errors thus depend on the physical location of dispersed photon events on the HRC-S, and 
so on both precise nominal pointing direction and the deliberate Lissajous figure "dither" 
about this point § . An updated version of the degap polynomial coefficients based on em- 
pirical wavelength corrections from multiple HRC-S/LETG observations of Capella has been 
released since CALDB 3.2.15. This update allows the correction of the nonlinearities in the 
HRC-S dispersion relation improving the 90% statistical significance uncertainties across the 
detector (i.e. 1.6a) from 22 mA to 16 mA. Since we have co-added three different obser- 
vations, and then the resulting positive and negative orders, we should consider a larger 
uncertainty. By propagating four times, in quadrature, the quoted deviations, we obtain 
90% uncertainties of 32 mA. However, since the updated wavelength corrections have been 
applied to the rest-frame position of the strongest soft X-ray metal electronic transitions and 
the interpolation of the correction may not be strictly valid to the position of blueshifted 
and redshifted lines, we conservatively adopt a 90% wavalength uncertainty 40 mA, 80% of 
the LETG Line Spread Function (LSF) FWHM. 

In the following we report all uncertainties at the la (68%) statistical significance level, 
with the exception of X-ray wavelength and redshifts, for which we use 90% statistical 
confidences (i.e. 1.6cr). 



3.2. Data Analysis 

We used the fitting package Sherpa, in CIAO, to fit our Chandra-LETG spectrum of 
1ES 1553+113. We ignored data outside the 10 — 60 A spectral interval, where the (S/N) 
was either too low (A > 60 A and A < 5 A) and/or the spectral resolving power too low (at 
A < 10 A, R< 200, implying Av > 1500 km s^ 1 ). 

We first fitted the spectrum with a model consisting of a single power law absorbed by 
neutral gas. This yielded a best-fitting power-law energy spectral index and cold-absorber 
column density of T = 2.66 ±0.01 and N H = (3.71 ±0.05) x 10 20 cm" 2 (fully consistent with 
the Galactic value along the line of sight to 1ES 1553+113: Ng ai = 3.65 x 10 20 cm" 2 (Kalberla 
et al. 2005). The integrated reduced x 2 was acceptable, xl{dof) = 0.90(3995). However, 
a visual inspection of the residuals showed both broad deviations around the instrumental 
edges of OI and CI and the HRC-S plate gaps (A ~ 40 — 54 A), as well as narrow deficit of 



See http://cxc.harvard.edu/cal/Letg/Hrc_disp/degap.html 



See, e.g. http://cxc.harvard.edu/cal/Letg/Corrlam/ 
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counts at specific wavalengths. Figure 1, top panel, shows these residuals, smoothed through 
a convolution with the LETG LSF, in re-normalized (after smoothing) standard deviations. 

At first, we excluded the possible absorption lines (which do not affect the continuum 
fitting over the integrated 10-60 A band) and focused on curing the residual calibration 
uncertainties at the position of the instrumental 01 and CI edges and the HRC-S plate gaps, 
with the aim of reaching an optimal model of the continuum that would allow us to properly 
assess the statistical significance of true physical spectral features possibly present in the 
spectrum. We then added to the continuum model, (a) an edge with negative optical depth, 
to cure the sharp edge-like emission feature present in the residuals at A ~ 43.8 A and due to 
either an over-estimated modeling of the CI instrumental edge in the LETG effective area, 
or a non-accurate modeling of the CI absorption edge in the wabs model of XSPEC, or both 
(best-fit r = —0.04); (b) two broad emission Gaussians at A = 41.61 A and A = 50.10 A 
(best-fitting o = 0.8 A for both) to model residuals broad wiggles at the wavelength of the 
CI edge, and the dithering-shaped HRC-S plate gaps ; (c) Four notch XSPEC component 
models: 3 negative (at A = 43.99, 48.34 and 52.14 A, with total widths of 0.4, 0.2 and 1.1 A, 
respectively), to model residual narrower (but still well resolved) deviations at the CI edge 
and plate-gaps wavlengths Fig. 1, top panel), and one positive with central A = 20.14 A 
and total width of 5.5 A, to model a systematic excess over the 18-23 A spectral region, just 
leftward the instrumental OI edge. 

We re-fitted the data, and obtained an improvement of A% 2 = 107, for 20 additional 
degrees of freedom, with almost unchanged best-fitting parameters: T = 2.62±0.01 and N# = 
(3.60 ± 0.05) x 10 20 cm" 2 . We measure a (0.1 - 2) keV flux of F .i- 2 = (2.45 ± 0.03) x 10~ n 
erg s _1 cm -2 (1.2 mCrab). The LSF-smoothed residuals now appear flat over the entire 
10-60 A band (Fig. 1, bottom panel). 

A visual inspection of the LSF-smoothed residuals in Fig. 1 (bottom panel), shows 
clear narrow (i.e. unresolved: AA<0.05 A) deficit of counts involving 3-6 consecutive bins, 
especially in the spectral region where redshifted C absorption is expected. This is expected: 
indeed, the sensitivity of dispersed spectrometers with constant AA and S/N per resolution 
element, over a given band-pass, to a given ion column density, increases with wavelengths. 
This is because, for a given column density, the equivalent width (EW) of a given absorption 
line scales with the square of the rest-frame position of the line, and linearly with the redshift. 



'During every pointing the Chandra satellite dithers around the aim-point sky coordinates producing a 
Lissajous-figure in detector coordinates. This implies that plate-gaps are not blind to photons at all times 
during the observation, and their dispersion position (for dispersed spectra) can be reconstructed a posteriori 
by applying the aspect-solution. However, the exact shape of the effective area at the plate-gap wavelengths 
depends on the accuracy of the aspect-reconstruction and the zeroth-order centroid. 
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For example, for equal ion column densities, the EWs of carbon absorption lines are about 
four times larger than that of corresponding oxygen lines. The Chandra - LETG effective 
area is similar in the redshifted O and C regions, therefore it is about 4 times easier detecting 
C than O with the Chandra-LETG (for equal C and O column densities). 

To identify candidate unresolved absorption lines we wrote a routine that scans the 
LSF-smoothed residuals (Fig. 1, bottom panel) over AA = 125 mA wide contiguous bins 
(2.5 LETG resolution elements), to search for negative deficits of counts involving between 
3-7 consecutive bins (i.e. 0.75-1.75 LETG LSF FWHM at our oversampled resolution), and 
line-shaped, i.e. symmetrically distributed around a central negative peak (at least one left 
and one right bin with values higher than the peak). We do not consider resolved lines here 
(i.e. consecutive deficits of counts spreading over more than ~ 7 bins), since we assume that 
any ISM or IGM absorption lines must have intrinisc FWHM <C 800 km s -1 . By definition, 
LSF-smoothing of the residuals tends to smooth flickering in adjacent sub-resolution bins. 
Thus, on the one hand it highlights real absorption (or emission) lines, but on the other hand 
it also increases the number of false absorption (or emission) line-like features with typical 
LSF width (see, e.g., Fig. 1 bottom panel). We can estimate the expected number of LSF- 
shaped emission or absorption features exceeding a bin-by-bin integrated (in quadrature) 
significance of 3a, by using a simple Gaussian argument. At our 4x LSF over-sampled 
binning, we have a total of (60-10)/0.0125 = 4000 bins. Therefore the Gaussian-predicted 
number of positive or negative LSF-shaped fluctuations at > 3a is (1 - 0.997) x4000 = 12. 
So, on average, we expect a total of 6 emission and 6 absorption LSF-shaped fluctuations. 
Indeed, the blind run of the scanning routine over the 10-60 A LSF-smoothed residuals 
(over-sampled at a resolution of 0.25 x the LETG LSF FWHM), finds exactly 6 LSF-shaped 
emission features exceeding the bin-by-bin integrated (in quadrature) statistical significance 
of 3a, according to expectations given the blazar nature of the target and the consequent 
absence of intrinsic emission in the optical or infrared spectra of 1ES 1553+113. On the 
contrary, 50 of such features are found in absorption, clearly suggesting that real absorption 
lines are imprinted on the 10-60 A LETG continuum spectrum of 1ES 1553+113. Given the 
way the scanning routine works, moving on contiguous 2.5 LETG-LSF wide bins, some of 
the 50 absorption line candidates found by the routine are simply the same strong feature, 
detected twice over adjacent scanned bins. This leaves about 40 absorption-line candidates, 
still many more than predicted by Gaussian statistics and found in emission. 

To confirm the real absorption-line nature of these candidates (at least down to the 
sensitivity of the current LETG spectrum), we therefore used unresolved (line FWHM frozen 
at 5 mA) absorption Gaussians to model the non-smoothed data. We added these Gaussians 
one by one and refitted the data, by allowing the line position to vary within 0.5 A around 
the scan-routine determined position of the centroid. This procedure yielded the results 
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summarized in Table 1, columns 1 and 2. Only 12 of the possible ~ 40 absorption candidates 
were actually confirmed by the fitting procedure at > 90% (1.6a) statistical significance. 
The single-line (prior to secure identification) statistical significances of 11 of these lines 
(computed as the ratio between the best-fitting line normalization and its la uncertainty) 
are in the range (2.2 — 4.1)a: three lines at > 4a, three at 3.6 < a < 4, five at 2.2 < a < 
2.8. (Table 1, column 2 and Fig. 1 bottom panel). An additional line-like feature is only 
hinted by the data and can be fitted with a negative Gaussian with single-line statistical 
significance of 1.7a. In the following we consider the six lines with prior-to-identification 
statistical significance > 3.6a (0.04 and 0.08 expected by chance up to z = 0.4 in the 
oxygen and carbon bands, respectively; see §4.2) as real absorption lines, which need to 
be identified. The remaining six absorption features are considered as tentative detections 
that need confirmation, and we consider them here only in association with either higher 
significance LETG absorbers or COS signposts, or both. In Figure 2 we show the best 
fitting models of the 12 absorption lines superposed on the non-LSF-smoothed residuals (in 
standard deviations) of the data to the best-fitting continuum model. Best fitting Gaussians 
models are color-grouped into 6 different systems, according to their possible identifications 
(§4-1). 



4. Discussion 



The LETG spectrum of 1ES 1553+113, reveals the presence of unresolved line-like 
absorption in both the z < 0.4 oxygen Ka and carbon Ka regions. These absorption lines 
could be the imprints of either Galactic ISM and/or CGM, or of warm-hot intervening 
IGM, or both. To identify these systems, we consider all the strongest ground-state metal 
transitions (Ka) from all ions of the two most abundant elements in gas with solar-like 
composition (O and C). We try to match the possible identifications, with the available 
redshift priors: (a) z = 0, rest-frame absorption by the ISM of our Galaxy, or CGM (the 
LETG resolution does not allow us to kinematically distinguish between the two; e.g. Gupta 
et al. 2012 and references therein), and (b) the redshifts of the possible BLAs (or NLAs) 
and/or OVI absorbers identified in the COS spectrum of 1ES 1553+113 (Danforth et al. 
2010). 



4.1. Absorption System Identifications 



The longest wavelength HI Lya securely identified in the COS spectrum of 1ES 1553+113, 
is at A ~ 1695 A, z — 0.395. At longer wavelengths, and up to the COS G160M band-limit 
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Table 1: Absorption Lines in the LETG (and COS) spectrum(a) of 1ES 1553+113 



"Wavelength 


b Significance Redshift 


Id 


COS HI, Metal I 


in A 


in a 








Galactic Absorbers 


23.34 ±0.04 


2.5 


±0.002 


OI Molecular 


N/A 


23.49 ±0.04 


4.1 


-0.001 ±0.002 


OI Ka Atomic 


N/A 


43.01 ±0.04 


2.8 


±0.001 


CII Ka Atomic 


N/A 


Intervening Absorbers 


High-Confidence X-ray System Detections and COS Associations 


44.27 ±0.04 


4.1 


0.312 ±0.001 


CVI Ka 


BLA-OVI z FUV = 0.3113 


52.82 ±0.04 


4.1 


0.312 ±0.001 


CV Ka 


BLA-OVI z FUV = 0.3113 


45.58 ±0.04 


3.8 


0.132 ±0.001 


CV Ka 


BLA z FUV = 0.1333 


38.25 ± 0.04 


2.7 


0.134 ±0.002 


CVI Ka 


BLA z FUV = 0.1333 


Tentative X-ray Detections and COS Associations 


49.80 ±0.04 


3.9 


0.237 ±0.001 


CV Ka 


NLA zfuv = 0.2356 


[ 26.65 ±0.04 


1.7 


0.234 ±0.002 


OVII Ka 


NLA zfuv = 0.2356 ] 


47.68 ±0.04 


3.6 


0.184 ±0.001 


CV Ka 


Triple-OVI z FUV = 0.186 - 0.190 


47.94 ±0.04 


2.2 


0.191 ±0.001 


CV Ka 


Triple-OVI z FUV = 0.186 - 0.190 


22.48 ±0.04 


2.3 


0.041 ± 0.002 


OVII Ka 


BLA zpuv = 0.0428 



a We use 40 mA as la systematic uncertainities for the HRC-LETG co-added negative and positive order 
wavelength scale. This is (4/5) of the LETG Line Spread Function FWHM, and it is due to the non-linearity 
in the HRC-LETG effective dispersion relation. 

b Single-line statistical significances, obtained by ratioing the best-fitting line normalization and its la error. 



13 



(A ~ 1800 A, or z ~ 0.48 in HI Lyot), the absorption line density drops dramatically, indi- 
cating z ~ 0.4 as a likely redshift for the target (Danforth et al. 2010). We therefore assume 
z = 0.4 as the redshift of 1ES 1553+113. 

The 12 absorption features identified as possible lines in the LETG spectrum of 1ES 1553+113, 
have wavelengths between A = 22.48 — 26.65 A (four lines in the O-region of the spectrum) 
and A = 38.25 — 52.82 A (8 lines in the C-region of the spectrum). The shortest wavelength 
Ka transition of O is that of its H-like ion (OVIII), at A = 18.97 A (Verner, Verner & 
Ferland, 1996), while the longest wavelength Ka transition of O is the inner-shell transition 
of its neutral ion (OI), at A = 23.52 A (E. Behar, private communication) . Similarly, for 
C, we have \(CVI La ) = 33.74 A (Verner, Verner & Ferland, 1996), and X(CI Ka ) = 43.37 
and 43.27 A (an LETG-resolved doublet; E. Behar, private communication), as the shortest 
and longest wavelength transitions, respectively. Intermediate ionization ions have their Ka 
transition wavelengths between these extremes . So, the four LSF-shaped features in the 
22-27 A range (Figure 1, top panel), if true absorption lines, can only be identified with 
OVIII in the ranges 0.185 < z < 0.4 or with OI-OVII at Otz < 0.4, with the exact lower 
limit depending on the ion considered. Analogously, the eight line-like features at A > 38.25 
A can only be identified with CVI in the range 0.134 < z < 0.4 or CI-CV at 0<z < 0.4, 
while the feature at A38.25 A, can only be due to redshifted CVI. 

We have checked these ions one by one, and propose possible identifications based on 
HI and/or OVI COS associations. We only check here for possible X-ray C and/or O coun- 
terparts to the FUV signposts. These are by far the two most abundant metals in gas with 
solar-like composition (Asplund et al. 2009), with C being only a factor of ~ 2 less abundant 
than O, and the third most abundant metals (N and Ne) being ~ 4x less abundant than C. 
Similarly, when in doubt between two or more identifications involving different ions of C 
and/or O, we always favor the He-like metal identification, i.e. CV and/or OVII. These are 
the He-like ions of C and O, respectively, and as such are highly stable and by far the most 
abundant ions of their respective elements over a broad range of temperatures. The oscillator 
strengths of their Ka transitions are also the strongest of all the other Ka transitions from 
lower- or higher-ionization ions. 



8 For most of the Ka inner-shell transitions of C and O (with the exception of OI and CII, for which 
empirical observational values are available, e.g. Takei, Fujimoto & Mitsuda, 2002, Shulz et al. 2003, Juett 
et al. 2004) we use wavelengths and oscillator strengths computed through the HULLAC code (Klapisch et 
al. 2003). The uncertainties on the rest-frame positions of these transitions are of the order of 0.1-1% (i.e. 
A A ~ 20 — 200 mA (i.e. Behar & Netzer, 2002). However, several of these positions have been bench-marked 
with either laboratory or astrophysical experiments (e.g. Kaastra et al. 2011, and reference therein; Shulz 
et al. 2003, and references therein; Juett et al. 2004)., and have resulted to be accurate to A A ~ 2 — 10 mA. 
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4-1.1. z ~ Galactic Absorption 

The three lines at A = 23.34 A, A = 23.49 A and A = 43.01 A, are consistent with 
z ~ absorption by neutral (or only moderately ionized) species of C and O. In particular, 
we clearly (4.1<r) detect atomic 01 at A = 23.49 ± 0.04 A and also a marginal (2.5<r) excess 
of molecular 01 (possibly ISM oxygen in coumpound forms of oxide dust grains, e.g. Takei, 
Fujimoto & Mitsuda, 2002) at A = 23.34±0.04 A, with respect to the instrumental molecular 
OI. Atomic OI is not an instrumental feature (e.g. Takei, Fujimoto & Mistuda 2002). Finally 
we marginally (2.8a) detect the signature of CII (whose Ka transition is twice as strong as 
the single transitions of the resolved and undetected CI doublet) at A = 43.01 ± 0.04 A. 

We identify these three transitions as due to the cold ISM of our Galaxy (Table 1, 
Columns 3-4; Fig. 1 and 2, yellow lines), and we therefore consider their single-line statistical 
significance as their true (i.e. after a secure identifications) statistical significance. 

We also note that we do not detect any OVII Ka z ~ absorber along the line of sight 
to 1ES 1553+113, down to the 3<r sensitivity limit of our Chandra LETG spectrum, in this 
spectral region: EW(OVII(^ = 0))< 12 mA. 

Below we present the possible intervening system identifications in order of statistical 
significances. 

4.1.2. The z x = 0.312 ± 0.001 CV-CVI-BLA-OVI WHIM System 

The line at A = 44.27 A, is that detected at the highest statistical significance (4.1<r) in 
the LETG spectrum (together with the line at A = 52.82 A and the Galactic OI atomic line). 
The CI, CII and CIV Ka transitions, at this wavelength, would be redshifted at z x = 0.021, 
0.031 and 0.070. At none of these redshifts is there a metal or HI association in COS. At the 
redshift of a putative CIII, instead, zx = 0.051, there is a moderately strong NLA. However, 
no imprint of a CIIIA977 absorber is seen in the FUSE spectrum of 1ES 1553+113 (though 
its presence cannot be ruled out due to the presence of HI Ly/3 and OI absorption/airglow 
in FUSE), and no other metal lines are present in the FUSE or COS spectra at this redshift. 
CV and CVI identifications are more promising. 

The CV Ka transition would be redshifted at z x = 0.100 ± 0.001. At a close (but 
> 3cr inconsistent) redshift (zpjjy = 0.10230) there is one of the strongest HI Lyct lines in 
COS, which is also seen in Ly/3 and is consistent with being thermally broadened (BLA) at 
T~ 10 5 K. Moreover, at A = 45.58 A there is another line in the LETG spectrum (3.8<r), 
which, if identified with CIV, would fall at zx = 0.100 + 0.001. The possible temperature of 
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the HI absorber at zfuv = 0.10230 is consistent with the presence of CIV and CV, however 
no CIV doublet (A = 1548.195,1550.77) is seen in the COS spectrum at z = 0.10230 or 
z = 0.100. We conclude (both for the lack of CIV in COS and the redshift inconsistency 
with the BLA) that the zx = 0.100 CV/CIV identification with a thick and moderately cool 
WHIM filament traced by the COS BLA at z FUV = 0.10230, is unlikely. 

A more likely identification is with a CVI Ka transition at A = 44.27 A, which would 
have a redshift of zx = 0.312. At A = 52.82 A there is the other highest statistical significance 
line of the LETG spectrum (4.1<r), which, if identified with CV Ka, would fall at zx = 
0.312 ± 0.001, exactly the redshift of the CVI identification. No HI or OVI identification is 
reported at a consistent redshift in the published COS spectrum. However, with additional 
data, the currently available COS spectrum of 1ES 1553+113 has a S/N a factor ~ 1.5 — 2x 
higher than the published spectrum. We therefore downloaded from the HST-MAST archive 
all available HST-COS spectra of 1ES 1553+113 taken with the G130M and G160M gratings, 
and re-analyzed the two portions of the spectrum at 1350 — 1365 A and 1590 — 1600 A, regions 
where the OVI(ls 2 2s ->■ ls 2 2p) doublet (i.e. OVIi(A1031.9), oscillator strength / = 0.133, 
and OVI 2 (A1037.6), / = 0.066) and the HI Lja at z ~ 0.308 - 0.316 are supposed to lie. 

We identify a strong HI BLA at z FUV = 0.31124 ± 0.00003 (8.1a; Fig. 3 top panel), 
with the following best-fitting parameters: A = 1594.03 ± 0.03 A, b = 58tg km s~ x and 
EW= 138 ± 17 mA. The strongest line of the OVI doublet (OVL) is also detected (3.9<r; 
Fig 3, central panel), with best fitting parameters: A = 1353.13 ± 0.04 A (i.e. zovi — 
0.31130 ± 0.0004), b = 38 ± 16 km s -1 , and EW= 27.5 ± 7.0 mA. Finally, the weakest line 
of the OVI doublet (OVI 2 ) is not detected in the data (Fig. 3, bottom panel), but the COS 
3o- EW upper limit at the wavelengths of this line, EW< 21 mA, is consistent with the 
maximum possible EW expected for the OVI 2 (given the measured EW(OVIi) and the fact 
that this line is not blank). We note that the HI/O VI Dopller parameter ratio is consistent 
with almost no turbulence contribution to the line breadth: b ~ 60 km s^ 1 for H, would 
correspond to logT= 5.32, at which temperature the OVI should be thermally broadened to 
b ~ 23 km s -1 , which is within la from the observed fe(OVIi). 

We conclude that the two highest significance intervening LETG lines are CV and CVI 
absorbers at zx = 0.312 (Table 1; Fig. 1 red lines), possibly imprinted by an intervening 
WHIM filament, as further suggested by the a-posteriori discover of an associated strong BLA 
and OVI counterpart in the COS spectrum, with the right temperature (i.e., line Doppler 
parameters) and strength. 
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4.1.3. The (z x ) = 0.133 ± 0.002 CV-CVI-BLA WHIM System 

The line at A = 45.58 A (3.8<r), cannot be associated with any metal or HI system in 
COS, if identified as CII, CIII or CVI. The CIV identification has already been ruled out in 
§4.1.2, because no corresponding CIV transition is detected in COS at z = 0.100 ± 0.002. 
At the redshift of a putative CI identification (z x = 0.051), instead, there is a NLA in COS 
(zfuv = 0.05094). However, CI must be associated with almost neutral gas, with a large 
HI/H fraction. The HI column derived from the COS data is only N HI = 2 x 10 13 cm~ 2 , 
implying a total baryon column of at most N# ~ 10 15 cm -2 . Galactic CI is not detected in 
these LETG data, and our line of sight to 1ES1553+113 crosses a Galactic column density 
Nfl- = 3.7 x 10 20 cm" 2 . Finally, there are 10 CI transitions accessible in the COS bandpass, 
and none is detected at z = 0.051. We therefore rule out the CI identification for this line. 

The only possibility left is CV at zx = 0.132 ± 0.001. The marginal detection at 
A = 38.25 A (2.7a) would have a consistent redshift of zx = 0.134 ± 0.002, if identified 
with CVI Ka. At a redshift consistent with both X-ray lines (zpuv — 0.13324) there is a 
weak (Nfl-/ = 10 13 - 15 cm~ 2 ) HI Lja in COS, with a doppler parameter of b = 46 ± 8 km s -1 . 
Correcting for the empirically derived b/b t h ratio of ~ 1.2 (Danforth, Stocke & Shull, 2010), 
gives a thermal width of b t h = 38^7 6 km s -1 , implying T = (0.6 — 1.7) x 10 5 K, consistent with 
either shock- heated WHIM or low-density purely photoionized IGM. However, we note that, 
given the relatively low-temperature derived from the HI Doppler parameter, the column 
density of this HI component is probably too low to produce all the CV and CVI detected in 
the LETG, even assuming Solar abundances. Therefore, if our identification is correct, the 
CV-CVI system cannot be entirely and directly associated with the HI Lja in COS. Possibly 
some of the detected CV and CVI must be associated with an even broader and shallower 
HI not detectable at the sensitivity of the current COS spectrum. 

We conclude that CV/CVI Ka at an average redshift of (zx) = 0.133±0.002. is the most 
likely identification for these two LETG lines (Table 1; Fig. 1, green line), possibly produced 
in a multi-phase WHIM filament together with the spatially associated BLA absorber. 

4.1.4. The z x = 0.237 ± 0.001 Multiphase WHIM System 

The line at A = 49.80 A is the second highest statistical significance line detected with 
the LETG (3.9a). This line can only be associated with a COS absorber if identified as a 
CV Ka at zx = 0.237 ± 0.001). Indeed, at a consistent redshift of zpuv — 0.23559, there 
is a weak NLA in COS. The HI line is too weak and much too narrow to be considered the 
actual counterpart of a putative CV absorber, so, if the LETG identification is correct, CV 
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and NLA would be the tracers of a spatially co- located multi-phase WHIM filament. 

We also note that, the lowest significance (1.7a) line-like absorption feature identified 
by our scanning routine and susbequently confirmed by a Gaussian fit in the LETG, lie at 
A = 26.95 A. If real, this feature would be identifiable with OVII Ka at z x = 0.234 ± 0.002, 
consistent with both the redshift of the putative CV absorber and that of the NLA in COS. 

We therefore tentatively identify the line at A = 49.80 A (3.9a) with a highly ionized 
shock-heated WHIM metal system at z x = 0.237 ± 0.001, with no detectable associated 
BLA at the sensitivity limit of the COS spectrum and embedded in a multi-phase system 
with a weakly photoionized portion producing the NLA at zfuv — 0.23559, possibly at the 
interface of the shock- heated part of the main-body filament (Table 1; Fig. 1, blue lines). 
Such multi-phase systems, with metals tracing hot collisionally ionized gas and detectable 
HI tracing a cooler photoionized portion of the same system, have been observed in NLA and 
OVI in the FUV (Danforth & Shull 2008), but it would be the first time they are observed 
in hotter gas traced by CV. Confirming the presence of the higher-ionization OVII absorber, 
together with CV and NLA, would probe an even more extreme regime and would therefore 
be extremely important. 



4.1.5. The z x = 0.184 - 0.191 Double-CV 

The line at A = 47.68 A (3.6a) can only be identified as CV Ka. No other ion of C 
can be associated with any metal and/or HI absorber in COS. Instead, the redshifts of the 
CV Ka identification, z x = 0.184 ± 0.001 is close, but only marginally consistent, with the 
short-wavelength extreme of a narrow redshift interval, z FUV = 0.18640 — 0.18989, where a 
rare triple-HI system is found in COS. No OVI is detected at zfuv = 0.18640, while both 
the zfuv = 0.18773 and the zfuv = 0.18989 absorbers have associated OVI absorption, 
with the central system (the one with higher column density) showing also NV, CIII and 
Sill absorption (Danforth et al. 2010). 

In the LETG we also marginally detect (2.2a) a line-like absorption feature at A = 47.94 
A. Also for this line, if real, the only possible identification, in association with a COS HI 
and/or metal system, is CV Ka at z x = 0.191 ± 0.001, consistent with the long- wavelength 
extreme of the COS triple-HI-metal system. 

We thus speculate that these two LETG lines may be the high-ionization CV coun- 
terparts (Table 1; Fig. 1, cyan lines) of the interesting, likely multi-phase, triple-HI-metal 
system detected in COS, possibly witnessing the on-going interplay between nearby galaxy 
outflows and a WHIM filament. 
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4.1.6. The z x = 0.041 ± 0.02 WHIM System 

Finally, the marginal detection at A = 22.47 A (2.3er) can only be redshifted OV or 
higher. The OV Ka transition would be redshifted to zx = 0.005 ± 0.002. There is a 
NLA present in COS at zpuv — 0.00717, i.e. marginally consistent with the X-ray redshift. 
However, OV is a tracer of photoionized gas that shows fractional abundances similar to OVI. 
No OVI absorption is seen at this redshift either in the FUSE spectrum of 1ES 1553+113 or 
in its LETG spectrum. The OVI Ka transition would be redshifted to z x = 0.020 ± 0.002. 
No HI Lya or OVI(ls 2 2s — > ls 2 2p) is seen at a consistent redshifts in COS or FUSE. 
Finally, the OVII and OVIII Ka transitions, would be redshifted to z x = 0.041 ± 0.002 
and Zx = 0.185 ± 0.002. These two redshifts are both consistent with either the broadest 
(hh — 55 — 80 km s _1 , corresponding to T= 10 5 ' 3 — 10 5 ' 7 K) and strongest (EW(BLA) 
= 135 ± 14 mA) BLA detected in the published COS spectrum, at zpuv = 0.04281, or 
the redshift range of the unusual triple-HI-metal system at zfuv = 0.18640 — 0.18989, 
respectively. However, while a possible double-CV counterpart to the triple-HI-metal system 
seen in COS is possibly detected in the LETG (§4.1.5), no OVII counterpart is detected. So, 
it appears unlikely that the line at A = 22.47 A is OVIII associated with the triple-HI-metal 
system. 

For this marginal LETG detection we therefore favor the z x = 0.041 OVII-BLA associ- 
ation (Table 1, Columns 3-5; Fig. 1 magenta line), and tentatively identify this association 
with an intervening WHIM system at the COS-detected BLA signpost (though an OVIII- 
Triple-HI association cannot be ruled out). 

We stress that this is the weakest system for which we propose an identification in 
association with a COS system (only one LETG line, at 2.3er). However, given the nature of 
this association (OVII-BLA), in principle probing the majority of the baryons in the WHIM, 
a confirmation of this putative OVII Ka detection at high statistical significance is highly 
desirable. 



4.2. On the Statistical Significance of the X-ray Absorption Lines 

The statistiscal significances of the X-ray absorption lines reported in Table 1 (second 
column), are single-line statistical significances. That is the ratio between the line flux, 
and its ler lower-bound uncertainty. For those lines not identifiable with z ~ transi- 
tions imprinted by Galactic or circumgalactic gas, these significances over-estimate the true 
absorption-line statistical significances, because they derive from a number of redshift trials. 
This number of redshift trials, depends on the wavelength of the transition, the spectrom- 
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eter resolution, and the available path length over which the blind search is performed: 
N= (A x Az)/AA spec , where A is the rest frame wavelength of the given transition, and 
AA spec the spectrometer LSF FWHM. In our case, Az = 0.4 and AA spec = 50 mA. So, the 
number of independent LETG resolution elements is 173 for the OVII Ka and 322 for the 
CV Ka. In such conditions, the true statistical significance (i.e. after correcting for the 
unknown line position) of a line blindly detected at a "single-line" statistical significance of 
4a is only 2.6a and 2.4a, for OVII Ka and CV Ka respectively. 

However, if an independent (and physically self-consistent) redshift prior is known (OVI 
or BLA in COS), then the line identification can be considered "secure" and therefore the 
single-line statistical significance provides a reliable estimate of the true statistical signifi- 
cance of the line. This is true for five out of the eight intervening lines: (a) the highest sta- 
tistical significance z x = 0.312 CV and CVI lines, at redshift consistent with that of a strong 
BLA-OVI system with the right temperature to produce CV, (b) the putative zx = 0.184 
and zx = 0.181 CV lines, at redshifts close to those of a clearly multi-phase triple-HI-metal 
system in COS, and (c) the speculative zx = 0.041 OVII absorber, at redshift consistent 
with that of a strong BLA whose thermal width is appropriate for the production of OVII. 

The remaining three lines are identified as CV Ka (and possibly OVII Ka) at zx = 0.237 
and CV-CVI Ka at (zx) = 0.133. Both lack a redshift prior (neither the NLA and the BLA 
at consistent COS redshifts, respectively, can be directly associated with the warm X-ray 
absorbers). However, for the (zx) = 0.133 system, the presence of two lines poses a bound on 
the number of redshift trials for the line detected at lower statistical significance and makes 
its measured single-line statistical significances close to the true statistical significances. The 
CV and CVI Ka at (z x ) = 0.133 are detected at single-line significances of 3.8a and 2.7a, 
corresponding to probabilities of chance detections of 1.5xl0~ 4 and 7 x 10" 3 , respectively. 
Allowing for the number of redshift trials, the probability of chance detection of the strongest 
of these two lines rises to 0.024, a Gaussian-equivalent significance of 2.3a. However, con- 
sidering this as a prior for the second line, the probability of having by chance both lines 
(independently) is the product of the two probabilities: P= 0.024 x 0.007 = 1.7 x 10 -4 , 
corresponding to a Gaussian-equivalent significance of 3.8<7. Finally, for the CV Ka line at 
z x = 0.237, the true statistical significance decreases from 3.9a (single-line significance) to 
2.4a, after allowing for the number of redshift trials (considering also the marginal OVII 
feature, the significance of the putative CV-OVII system would rise to 3.1a). 



-20 - 



5. Conclusions 

We report on the detection of six absorption lines and five additional features associated 
with either other X-ray lines detected at higher significance and/or Far-Ultraviolet (FUV) 
signposts, in the 500 ks Chandra-LETG spectrum of the bright (F ai _ 2 = 1.2 mCrab) and 
distant (z > 0.4) blazar 1ES 1553+113. The six detections have single-line statistical signifi- 
cances between 3.6 — 4. la, while the five additional associations are detected with single-line 
statistical significances of 2. 2-2. 8a. 

We tentatively identify these absorption lines as belonging to six different systems: the 
Galaxy cold ISM, and five intervening warm-hot IGM filaments. 

Two out of the five intervening absorption systems are detected in two different ions at 
true X-ray statistical significances of 5.8a (zx = 0.312) and 3.8a ((zx) = 0.133), and both 
have a BLA (and OVI for the system at z x = 0.312) counterpart at consistent redshifts in 
the HST COS spectrum of 1ES 1553+113. For these two systems, the COS counterparts are 
also physically consistent with the ionization degree and temperature of the X-ray absorbers 
in a shock-heated gas scenario. However, for the (zx) = 0.133 system, the BLA column 
density is probably too low to produce all the detected CV and CVI absorption. 

The three remaining proposed identifications are X-ray single-ion tentative detections, 
all with a possible, physically consistent, COS counterpart: (a) the double CV Ka system 
at z x = 0.184 (3.6a) and z x = 0.191 (2.2a), both with NLA (and OVI for the putative CV 
at z x = 0.191) counterparts, (b) the possible OVII Ka absorber at z x = 0.041 (2.3a) with a 
BLA counterpart, and (c) CV Ka absorber at at z x = 0.237 with an NLA counterpart (and 
possibly associated OVII absorption). 

Here we summarize our identifications: 

• Galactic ISM: We detect atomic OI (and possibly an excess of ISM oxygen in compound 
forms, over the instrumental molecular OI), as well as CII at z ~ and associate these 
absorbers with the cold ISM of our Galaxy along the line of sight to 1ES 1553+113. 

• Intervening System I: this is the highest significance LETG system, with two lines 
detected at single-line statistical significances of 4.1a each. We identify these two lines 
with CV and CVI absorption at z x = 0.312 + 0.001, possibly imprinted by an interven- 
ing WHIM filament, as further suggested by our a posteriori discover of an associated 
strong BLA and OVI counterpart in the COS spectrum, with the right temperature 
(i.e. line Doppler parameters) and strength. The combined (in quadrature) statistical 
significance of this system in the LETG is 5.8a. 
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• Intervening System II: we identify a CV-CVI-BLA system at (zx) = 0.133 ± 0.002, 
and associate this absorber with an intervening warm WHIM filament, where only 
part of the X-ray absorption is physically related to the detected BLA. The statistical 
significance of this system in the LETG (after properly accounting for the number of 
redshift trials) is 3.8<x. 

• Tentative Intervening System III: we identify 2 CV Ka absorbers at redshifts z x = 
0.184 ± 0.001 and zx = 0.191 ± 0.001, adjacent to the redshifts of the extremes of 
a narrow redshift interval containing a rare triple-HI-metal absorption system. We 
speculate that the X-ray CV absorbers trace the same multi-phase medium traced 
by the triple-HI-metal system in COS, possibly representing the result of the interplay 
between galaxy winds and an IGM filament. The two lines of CV are detected at single- 
line statistical significances of 3.6<x and 2.2<r (close to the true statistical significances 
because of the HI-OVI priors), respectively. 

• Tentative Intervening System IV: we identify CV (and possibly OVII) absorption at 
zx = 0.237 ±0.001, consistent with the redshift of a weak NLA. We speculate that the 
NLA be considered as an indirect (i.e., not physically associated) association tracing 
a multi-phase WHIM filament, with a hot, shock- heated part producing the X-ray 
absorption and a cooler photoionized part giving rise to the narrow HI absorption. 
The true statistical significance of this line in the LETG, allowing for the number 
of redshift trials, is 2.4<r (versus a single-line statistical significance of 3.9<r). If the 
association with the OVII is considered, the statistical significance of the CV-OVII 
system is 3.1a. 

• Tentative Intervening System V: we tentatively identify an OVII-BLA system at z x = 
0.041 ±0.002 and associate this absorber with an intervening hot WHIM filament. This 
is the weakest system for which we report a possible identification, with only one X-ray 
line at a single-line statistical significance of 2.3<r, but with a possible redshift prior 
identified in the broadest and strongest BLA reported in the published COS spectrum. 

This moderate quality 500 ks Chandra-LETG spectrum of the, by far, brightest steady 
soft X-ray target in the z > 0.4 sky, strongly hints at the existence of possible X-ray metal 
WHIM counterparts to genuine gaseous signposts, previously secured with the HST-COS 
(or viceversa), and so demonstrates the efficiency of this observational strategy. Moreover, 
these results suggest that the gas may be often found to be multiphase, possibly indicating 
that the IGM structure is richer than typically found in hydrodynamical simulations. Longer 
integrations with both the high resolution soft X-ray spectrometers available (the Chandra- 
LETG and the XMM-./Vewfon-RGSs) are needed to confirm at high statistical significances 
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and with two indipendent measurements these tentative identifications and associations (es- 
pecially for the three lowest statistical-significance associations), and to conduct a sensitive 
exploration of even hotter WHIM systems along this line of sight. 

6. Acknowledgements 

We thank Ehud Behar for providing wavelength, oscillator strengths and transition 
probabilities of inner shell transitions of C, and Ne. FN acknowledges support from 
CXO grant G01-12174X, and from INAF-ASI ADAE grant 1.05.04.13.15. MS and CD 
acknowledge support from COS grant NNX08AC14G to the University of Colorado, from 
NASA. XB acknowledges support by the Spanish Ministry of Economy and Competitiveness 
through grant AYA2010-21490-C02-01. 

REFERENCES 

Asplund, M., Grevesse, N., Sauval, A.J., Scott, P., 2009, ARA&A, 47, 481 
Bennett, C.L. et al, 2003, ApJS, 148, 1 

Bertone S., Schaye J., Dalla Vecchia C, Booth C. M., Theuns T., Wiersma R. P. C, 2010, 
MNRAS, 407, 544 

Branchini, E., et al., 2009, ApJ, 697, 328 
Brinkman, A.C., et al., 2000, ApJ, 530, L111B 

Buote, DA., Zappacosta, L., Fang, T., Humphrey, P. J., Gastaldello, F., Tagliaferri, C, 
2009, ApJ, 695, 1351 

Cen, R. & Ostriker, J.P, 2006, ApJ, 650, 560 
Danforth, C.W. et al., 2011, ApJ, 743, 18 
Danforth, C.W. et al., 2010, ApJ, 720, 976 

Danforth, C.W., Stocke, J.T. & Shull, J.M., 2010, ApJ, 720, 976 

Danforth, C.W. & Shull, J.M., 2008, ApJ, 679, 194 

Fang, T., Canizares, C.R., & Yao, Y., 2007, ApJ, 670, 992 

Gupta, A. et al., 2012, ApJ, 756, L8 

den Herder, JW. et al., 2001, A&A 365, L7 

Juett, A.M., Schulz, N.S. & Chakrabarty, D., 2004, ApJ, 612, 308 
Kaastra, j.S. et al., 2011, A&A, 534, 37 

Kaastra, J.S., Werner, N., den Herder, J.W.A., Paerels, F.B.S., de Plaa, J., Rasmussen, 
A.P., de Vries, CP, 2006, ApJ, 652, 189: K06 
Kalberla, PMW. et al., 2005, A&A, 440, 775 



-23 - 



Keshet, U., Kushnir, D., Loeb, A., Waman, E., 2012. larXiv: 1210. 15741 

Kirkman, D. et al, 2003, ApJS, 149, 1 

Mathur, S., Weinberg, D.H. & Chen, X., 2003, ApJ, 582, 82 

Nicastro, F. et al, 2010, ApJ, 715, 854 

Nicastro, F., Mathur, S. & Elvis, M., 2008, Science, 319, 55 

Nicastro, F. et al, 2005a, Nature, 433, 495 

Nicastro, F. et al, 2005b, ApJ, 629, 700 

Nicastro, F. et al, 2002, ApJ, 573, 157 

Murray, S.S. & Chappell, J.H., 1995, SPIE, 597, 274 

Rasmussen, A. P., Kahn, S.M., Paerels, F., den Herder, J.W., Kaastra, J., de Vries, C, 2007, 

ApJ, 656, 129: R07 

Rauch, M., ARA&A, 36, 267 

Schulz, N.S., Juett, A., Chakrabarty, D., Canizares, C.R., 2003, AN, 324, 166 

Shull, J.M., Smith, B.D., Danforth, C.W., 2012, ApJ, 759, 23 

Smith, B. D, Hallman E. J., Shull J. M., O'Shea B. W., 2011, ApJ, 731, 6 

Spergel, D.N. et al., 2007, ApJS, 170, 377 

Takei, Y., Fujimoto, R. & Mitsuda, K., 2002, ApJ, 581, 307 

Tornatore L., Borgani S., Viel M., Springel V., 2010, MNRAS, 402, 1911 

Treves, A., Falomo, R. & Uslenghi, M., 2007, A&A, 473, L17 

Verner, D.A., Verner, E.M. & Ferland, G.J., 1996, Atom.Ph, 4003 

Weinberg, D.H., 1997, ApJ, 490, 564 

Williams, R.J., Mathur, S., Nicastro, F., Elvis, M., 2006, ApJ, 624, L95 
Williams, R.J, Mulchaey, J.S., Kollmeier, J.A. & Cox, T.J., 2010, ApJ, 724, L25 
Yao, Y., Shull, J.M.„ Wang, Q.D., Cash, W., 2012, ApJ, 746, 166 
Yoshikawa, K. et al., 2003, PAS J, 55, 879 

Zappacosta, L., Nicastro, F., Krongold, Y., Maiolini, R., 2012, ApJ, 753, 137 
Zappacosta, L. et al., 2010, ApJ, 717, 74 



This preprint was prepared with the A AS IATgX macros v5.2. 



b 
G 

■ i— i 

Us 
P 

K 



-4 - 



20 



_i i i_ 



30 



40 



50 




-2 " 



-4 - 



20 



30 40 

Wavelength (in A) 



50 



Fig. 1. — LSF-smoothed residuals (in re- normalized, after smoothing, standard deviations) of the 20-55 A 
500 ks LETG spectrum of 1ES 1553+113, to: (a) a model consisting of a power-law absorbed by a neutral 
absorber (top panel), and (b) the same continuum model with the addition of 7 components (an edge with 
negative optical depth, two broad Gaussians, three negative and one positive notches, to cure the residual 
systematic calibration uncertainties visible in the top panel and due to the instrumental OI and CI edges 
and the dithering-shaped HRC-S-LETG effective area at the wavelengths of the HRC-S plate gaps (bottom 
panel). In the bottom panel, the positions of the 12 possible absorption lines are marked with color-coded 
segments (§4) and labeled with their single-line statistical significances. 
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Fig. 2. — Non LSF-smoothed residuals (in a) of four portions of the 500 ks LETG spectrum of 
1ES 1553+113, where Galactic OI and CII, as well as intervening OVII and CV, are tentatively detected, 
at single-line statistical significances in the 2.2 — 4.1c range. Best-fitting curves of different colors identify 
different systems: (1) Galactic OI and CII (yellow), (2) OVII at zx = 0.041 (magenta), (3) CV (and myabe 
OVII) at z x = 0.237 (blue), (4) CVI and CV at z x = 0.312 (red), (5) CV and CVI at (z x ) = 0.133 (green), 
and (6) CV at the COS triple- HI- metal redshifts (cyan). 
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Fig. 3. — Three portions of the COS spectrum of 1ES 1553+113, centered around the HI Lya (top panel), 
OVIi (central panel), and OVI2 (bottom panel), at Zfuv — 0.31124 ± 0.00003. In all panels, red curves are 
the best-fitting Gaussians for the HI Lya (top panel), OVIi (central panel) and OVI2 (bottom panel). COS 
data are binned at about 1/8 the COS resolution, i.e. A A ~ 9 mA. 



